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Q^ . Abstract. 

Broad band spectra of total and compact-scale radio emission from 
^^ • blazars, used in combination with kinematic information inferred from 

li I VLBI monitoring programs, can be applied successfully to better con- 

strain the models of compact jets. We discuss here a "hands-on" ap- 
lO ! proach for tying together the kinematic and spectral properties of radio 

emission from blazars. 
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O ' 1. Introduction 

(N 

^^ ' Analytical models of parsec-scale jets, although being contested ever stronger 

Q>^ , by numerical simulations, still provide a valuable means for describing obser- 

vational data, particularly when different aspects of the jet physics (such as 
Q^. its spectral and kinematic properties) are combined together in a single formu- 

lation, providing additional constraints and checks for the model parameters. 
In this contribution, we discuss how the well-established shock-in-jet model 
("shock model" hereafter; see Marscher 1990, Marscher, Gear, & Travis 1991, 
Gj ■ for detailed discussions of the model) can be reinforced by inclusion of kinematic 

information available from VLBI observations. 
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H ' 2. Model quantities 



In its most common formulation, the shock model predicts changes of the turnover 
frequency, fm, and flux density, S^a, in the spectrum of radio emission associated 
with a shock. The jet is usually assumed to have a constant opening angle, (/>, 
so that the shock transverse dimension is proportional to the distance, r, at 
which the shock is located. Other model parameters are expressed as functions 
of r: the magnetic field B oc r~'^, Doppler factor S oc r , and number den- 
sity A^ oc r~" (for a power-law electron energy distribution, N{^)d'y oc 'y~^d'y). 
The shock emission is dominated subsequently by Compton, synchrotron and 
adiabatic losses. At each stage, the predicted quantities are described by the 
following proportionalities: 5m oc i/^ and fm oc r^, with p = p{a,b,s) and 
£ = e(a, b, s) (for complete evaluations of p and e, see Marscher 1990, Lobanov 
& Zensus 1999). Below, we describe how estimates of the power index b can 
be obtained from VLBI data, and give examples of applying this approach to 
studying compact variable radio sources. 



3. Observable quantities 

Single dish observations yield light curves S{t) at different frequencies, which 
can be used to determine the evolution of spectral turnover (5*111, ^'m), provided 
an adequate frequency coverage and time sampling. VLBI monitoring programs 
allow to measure relative proper motions, ^app(i)) and (in exceptional cases) 
also spectral changes of enhanced emission regions detected in parsec-scale jets. 
From the measured /i(t), the jet apparent speeds, (3a_pp{t), and Doppler factors, 
6{t), can be reconstructed, with necessary assumptions made about the jet kine- 
matics. In the simplest case, the jet Lorentz factor, 7, can be taken constant, 
and 6{t) is then described by changes of the jet viewing angle, 0{t). For more 
complicated cases, 7(t) = 7mm(0 = [1 + /3app(0]°'^ can be assumed, or even 
complete kinematic settings can be postulated (e.g. a helical trajectory, as has 
been done, for instance, by Roland et al. 1994). 

4. Relations between the jet spectrum and kinematics 

Once the form of 6{t) has been determined, we can evaluate b{t). Since variations 
of 6{t) are not necessarily monotonic, we resort to determining b{t) locally, so 
that 

^ ' log[r(t + dt)/r{t)] ^ ' 

We then select a timerange, (^1,^2)5 during which the changes of h{t) are small 
enough to approximate h{ti) « 6(^2) ~ h. Fitting the observed spectral turnover 
data, we can obtain the turnover points at the respective epochs, and evaluate 
the absolute location of the shock at the epoch ty. 

with At = t2 — ti, Tu = (t'm2/^mi)- Repeating this step as many times as 
necessary, we can reconstruct the entire kinematic evolution of the shock. The 
procedure can also be reversed: we can first fit the shock model to the spectral 
data, and determine values of h for different time periods. We then use equation 
(2) to calculate the respective locations of the shock at different epochs, and 
compare these locations with the locations and speeds inferred from VLBI data. 
If we fix the kinematic settings and assume that the shocked feature moves 
at a speed /3j along a helical path with amplitude, A{r)^ frequency, u), and parallel 
wavenumber, /c, we can reconstruct the time evolution of the shock location: 

t{r) = to + 1^ fc^^2(^) + [C2(r)/32 -u;2fc2Ci(r)]V2^'' ' ^^^ 

with Ci(r) = [yl^(r)]2 + 1 and C2(r) = Ci(r) + k'^A^{r). The form of A{r) may 
differ, depending on the choice of the jet geometry. We use A{r) = A(ro)r/(ao + 
r), corresponding to a jet with opening half-angle approaching arctan[^(ro)], 
for r ^ tq. The obtained t{r) can be then checked against 6(t) inferred from 
the shock model, or used for predicting the light curves directly. 
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Figure 1. Observed (left) and simulated (right) lightcurves of a 
short-timescale flare in 0235+164. In the observed lightcurves, under- 
lying emission has been subtracted. The model lightcurves reproduce 
the observed timelags, and represent only the first of the two events 
(J.D. 8901-8907) seen in the observed light curves. 



5. Examples 

We show here two examples of applying the method outlined above to radio 
observations of blazars. 



5.1. 0235+164 

A short-timescale flare in October 1992 was monitored with the VLA at 1.4, 
4.8, and 8.4 GHz (Kraus, Quirrenbach, Lobanov, et al., these proceedings) . 
The observed light curves show that the emission first peaks at 1.4 GHz, and 
later nearly simultaneously at 4.8 and 8.4 GHz. A cross correlation function 
analysis yields the respective time lags r^'g = 0.8 + 0.2 days, T^;f = 0.7+0.2 days, 
and Tg'l = — 0.2 + 0.2 days. The flare duration becomes progressively longer at 
higher frequencies, making this event rather peculiar. The modified lightcurves 
obtained after subtraction of the underlying emission are shown in the left panel 
of Figure [l|. 

We discuss elsewhere (Kraus et al. 1999) several possible schemes capable 
of explaining the observed peculiarities. One of our proposed schemes uses 
a precessing electron-positron beam (see Roland et al. 1994) with a period 
Pq = 200 days and precession angle fig = 5.7°. The kinematics of such a beam is 
then given by equation (3), with A(ro) = 0.1 pc and oq = A{ro)/ tanilo = Ipc. 
The resulting Doppler factors vary with time, reproducing the observed lags 
between the peaks in the lightcurves. For time lags to be present during a flare, 
the turnover frequency in the observer's frame should be within the range of 
observing frequencies. In 0235+164, we can satisfy this condition by postulating 
a homogeneous synchrotron spectrum with spectral index a = —0.5 and rest 
frame turnover frequency u'^ = 0.15 GHz. Additional spectral evolution may be 
required to remove the apparent discrepancy between the model and observed 
amplitudes and widths of the flares. 
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Figure 2. Kinematic and spectral evolution of an enhanced emission 
region C5 in the jet of 3C 345. Left panel shows the observed path of 
C5 in the plane of the sky; right panel presents the evolution of the 
spectral turnover. Solid line in the right plane shows a fit by the shock 
model. Dashed line shows how the fit changes if the spectral evolution 
of the shock is required to comply with the observed path of C5. 

5.2. 3C345 

We have studied (Lobanov & Zensus 1999) spectral changes in the core and 
several jet components in 3C 345, based on the data from a VLBI monitoring 
of the source. In the example shown in figure ^, we use the observed trajectory 
(left panel of fig. ||) of the jet component C5 to confront a fit by the shock 
model to the variations of Sm and v^ (right panel of fig. y). In the right panel, 
the solid line shows a fit by the shock model, without taking into account the 
observed kinematics of C5. When we require the shock model to reproduce b{t) 
needed to satisfy the observed path of C5, the fit becomes problematic, at later 
stages of the shock evolution. This indicates that, at distances > 1 mas, the 
shock may have dissipated, and other processes become main contributors to 
the emission from C5. Incidentally, the kinematic and emission properties of 
other jet features also show evidence for a change in the emission properties, at 
distances of 1-1.5 mas from the VLBI core of 3C345 (Lobanov &: Zensus). 
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